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The Unear 5 expansion is used to obtain corrections up to 0{5^) to the self-energy for a complex 
scalar field theory with a \{ip*ip)'^ interaction at high temperature and non-zero charge density. The 
calculation is done in the imaginary-time formalism via the Hamiltonian form of the path integral. 
Nonperturbative results are generated by a systematic order by order variational procedure and the 
dependence of the critical temperature on the chemical potential fi is obtained. 
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. . I. INTRODUCTION 

o : 

I It has been known for some time that conventional perturbation theory is inadequate for describing high-temperature 
^vq r field theory [Q , with the perturbation expansion breaking down at some order in the coupling constant in the parameter 
^ ' regime where the temperature is very much greater than the bare mass. There have been several attempts to 
circumvent this problem (for example Refs. |^,^ use resummed perturbative expansions and systematically include all 
■ relevant diagrams) . A temperature-dependent renormalization group approach has been used in Ref . jj] . The method 
^—i , we will employ in this paper is the linear S expansion (see, for example Q), which contains elements of both of the 
QQ ' above methods. It involves a mass shift which is determined order-by-order in the expansion by a non-perturbative 
I criterion, but prior to the (crucial) optimization stage it merely uses low-order perturbation theory with modified 
' propagators and vertices. 

^ ' Pinto and Ramos joj have successfully applied the method to finite-temperature symmetry breaking in a scalar field 
^ theory with a ip'^ interaction up to second order in S, which in particular involves evaluating the non-trivial "setting- 
sun" diagram. In the present paper we extend this approach to tackle the problem of a non-zero chemical potential; 
in so doing we consider a complex scalar field rather than a real scalar field. For — 0, the effect is simply to alter 
the symmetry factors of the one- and two- loop diagrams, but for fi ^ there is a shift in the energy component of the 
' Matsubara propagator, which makes the diagrams considerably more difficult to evaluate but presents no problem of 
1 principle. 

(-H ■ To our knowledge there have been only a few attempts at this problem. Benson et al. performed a one-loop 
. calculation, while Funakubo and Sakamoto used a renormalization group approach in the large- A'^ limit of an 
O I' 0{N) theory (in which the setting-sun diagram does not appear). In the standard lattice Monte-Carlo approach it 
, is impossible to incorporate a non-zero chemical potential, because the Euclidean action then becomes complex and 
i-C ' cannot be used as a real, positive statistical weight. 
^ ] We outline the linear S expansion technique below and explain how it can generate nonperturbative results, while 
in Sec. II we present briefly the construction of the theory in the imaginary time formalism and include the chemical 
potential associated with a conserved charge. We then go on to calculate all diagrams contributing to the self-energy, 
^ . up to second order in S. Sec. Ill includes results for the critical temperature and its behaviour as one allows fj, to 
■ ■ ' take non-zero values. The appendix contains a detailed calculation of the two-loop setting sun diagram in the high 
temperature expansion. 

A. The Linear Delta Expansion 

The linear 6 expansion (LDE) is a technique that allows the use of an analytic approach to probe the nonperturbative 
sector of the field theory to which it is applied. It has been employed with success in a wide variety of areas 
with convergence of the expansion rigorously demonstrated in some simple zero- and one-dimensional models (see, for 
example, [^). The method involves the introduction of an artificial expansion parameter, S, and the modification of 
the action of the theory under consideration via 

s ^ s' = {1 - s)So{{m}) + 5S, (1) 
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where {rji} is some set of variational parameters and 5*0 represents the action of some soluble theory. This trial action 
is not determined by the method; however, an appropriate choice is usually suggested by the form of the theory under 
consideration. The (5-modified action can then be used to evaluate some desired physical quantity as a power series in 
d (truncated at some finite order, N), with 6 then set equal to 1 at the end of the calculation. We label this quantity 
P/V) noting that it will, in general, have a residual dependence on the {rji}. The variational element is introduced by 
fixing these parameters according to some specified criterion. Wc shall use the principle of minimal sensitivity (PMS) 
p3 whereby the {r]i} are chosen at a stationary point of the quantity Pn, namely: 



dP, 



N 



dm 



= 0. (2) 

rii=rii 



It is this order by order fixing that allows for non-perturbative behaviour to emerge (the optimized variational pa- 
rameters becoming functions of the order of truncation, and so being more correctly labelled by {?7i(-^V)}), and can 
also provide for convergence of the expansion, two particularly desirable features. 



II. THE CHARGED SCALAR FIELD AT FINITE TEMPERATURE 



In Minkowski space, the Lagrangian density of a massive, complex scalar field ip{x) with a V{ip*ip) interaction term 
is given by 



(3) 



This Lagrangian density has a well-known global U{1) gauge symmetry which leads to a conserved charge Q and an 
associated chemical potential /u, so that the grand partition function is 



Z(/3,m) =Tre-^(^-'^'5). 



(4) 



Proceeding through the Hamiltonian form of the path integral for Z, in the imaginary-time formalism of finite- 
temperature field theory p2|, we arrive at 



where the Euclideanized action is (dropping the suffix for future notational convenience) 

S{j3, fi) = Sf(P, /i) + Si{f3, /i) + Sc{f3, m), 

with 



d"^ d 



Siip,^i) = / d\v(^V), 

JT 



-4( g7I-2/i^+V- + M' 



ip+ / Acv((^*^), 



(5) 

(6) 

(7) 
(8) 
(9) 



where t — it, x — (t, x) and Jrpd'^x — dr Jd'^x. Here Sc{P,tJ) represents those counterterms required to render the 
model finite. 

A sensible choice for our delta-modified action, in the presence of a {tf*LpY self-interaction, V{Lp*Lp) — {\/ A){if* LpY , 
would then be 



S\P, Ai) = S%{P, /i) + SfiP, Ai) + S'ciP, t^), 



with 



^f(/3,m) 



d'^x ip* 



d_ 



d'^x 



SX 



d'^x ip* 



-A' 



2/1 -- + + /i^ 







/ d^x 






JT 



(10) 

(11) 

(12) 
(13) 
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introducing the variational parameter fi, or equivalently 77, given by i]'^ = fl^ — w?. In subsequent sections we will 
evaluate the thermal mass up to 0(15^), using dimensional regularization to renormalize the theory. It is an important 
point, raised in ||^, that the arbitrary variational parameter 77 becomes a function of the bare parameters (which is 
how non-perturbative behaviour arises in the LDE), so we must eliminate any divergences before we apply the PMS 
optimization procedure. The renormalization procedure we use is identical to that of |^ and is based on |TT| ; we forgo 
any detailed discussion of the cancelling of temperature-dependent divergences and any systematic calculation of the 
coefficients of the counterterms , and C"^, the details being essentially identical to those discussed in 



A. Calculating the thermal mass 



We recall that the Feynman rules in frequency space in the presence of an overall charge are as follows: 

1. To every line of the diagram assign a factor Ai?(ia;„,k) = [—{iujn + /i)^ + lo\]^^ , where ojn = 27rn//3 and 

— \/(k^ + r2^), and an arrow in the direction of momentimr flow 

2. Assign a factor —5\ to each vertex 

3. Assign a factor Srf' to each insertion 

4. Integrate over every internal line with the measure Tj^nl d^^/{2TT)^ 

5. There must be conservation of charge at each vertex, i. e. the number of arrows entering a vertex must be the 
same as the number leaving 

We are now in a position to estimate corrections to the thermal mass (which we do up to 0{S'^)), defined by 

m^^^ = n^ +U{iLJn,p), (14) 

evaluated on-shell, i.e. iuJn = — /x, p = 0, where 11 is the thermal self-energy. This choice of four-momentum is 
discussed in the Appendix, where the setting-sun diagram is calculated, the only momentum-dependent contribution 
to n up to 0{S^). 



B. The thermal mass at first order 



1 2 
FIG. 1. Diagrams contributing at 0{5) 

To lowest order, the relevant contributions are 

= -Sf^'' + S\TY, / T^AF{tLOn,k). (15) 

The frequency sum in ( |l5|) can be treated in a particularly concise and efficient manner through the mixed represen- 
tation ||l| AF(T,k) of Af: 

Af{t, k) = T ^ e-'^-^ApiiLOr,, k) 

n 

= [(1 + )e-("''-'')" + n+e^'^^+^^^] , (16) 



where 



3 



is the Bose-Einstein distribution function in the presence of a non-zero chemical potential. The Matsubara propagator 
is recovered by Fourier transforming the mixed propagator with respect to the r variable: 



Jo 

The contribution of IIj in can then be calculated trivially to give 

Using dimensional regularization [Q one finds that the 0{6) contribution to the thermal mass is 



(18) 



1 r \ 



where M is a mass scale introduced by dimensional regularization and, in the notation of 

>T'l{y,r) ^ ^[h3{y,r) + {y, -r)] , 



hsiy.r) 



dx- 



r(3) Jo V^^T? 



exp{^yx^ + y^ — ry) — 1 



(19) 

(20) 

(21) 
(22) 



where x = (3k, y — /3r2 and r — ji/Q,. Using the minimal subtraction scheme (MS) we can eliminate the divergent 
term arising from the temperature-independent part of the self energy using the 0{5) mass counterterm 



We now analyze (|2^) in the limit y ^ (the high temperature expansion) keeping r fixed jl^ : 



h^(v,r) = -\/l — — — 



leading to the following expression for the renormalized thermal mass at first order: 



r^2 . 2 .Ar2 Ara r- 



12 ' 4tt " ^ ' ^ 167r2 
In this and all subsequent calculations we neglect terms of 0{1/T). 



In 



7- 2r^ 



(23) 



(24) 



(25) 



C. The thermal mass at second order 
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At C((5^) there are five diagrams (shown in Fig. 2) which provide corrections to the self-energy. We begin by 
evaluating the first diagram, which in the high temperature limit is 



nf 



1 



6' 



2 W 
167r2 



( 



(26) 



The next two diagrams are constructed using mass and vertex counterterms to give 0{S^) diagrams that eliminate 
the temperature-dependent divergences arising from . Explicitly, we obtain in the high temperature limit 



n 



(167r2)2e2 
202 



6'- 



A2 



1 



2(167r2)2 

, 5A^r>2 
2(167r2)2e2 

, 5A2f]2 



In 



V47rM2 J 



8n yT 

1 2 

+ 7 4 



167r2 



5A2 



rp2 

327r2e 1 12 



In 



V47rAf2 



4tt 
-I 2 



Vl - r2 + 



+ 7-1 



(327r2)2 

The momentum-independent two-loop diagram is given by 

„o A2rj2 1 „^ a2 1 rr2 to , m i 



1 + 



(167r2 



167r2 e 1 12 47r ^ 



Stt Vl - r2 



A2T3 



1 



967rf7 ^1 - r2 
47rT2- 



In 



V M2 



■7 



A2T2 
32^ + 
A2TO 



. A2 

(8^)2 
2 



^ Vi 

3 TT 



In 



r!2 



llVf2" 



167r2 



In 



/47rT2 



87r2 
TVl 



In 



/47rT2 

(vl^72~ 



7- 2r^ 



27r 



ill 

47r2 



7 — r 



In 



V Af2 



647r3 Vl-r2 



(52 



A202 
(167r2)2 



(27) 



(28) 



7 - 2r^ 



(29) 



Finally, we have the momentum-dependent setting sun diagram, which we evaluate with the Euclidean self-energy 
on shell (the details of the calculation are provided in Appendix A) : 



n = 



3X^n^ 1 

(327r2)2^ 
3A2r!2 



^9 3A^r2^ 1 _o 



A2p2 



(327r2)2 e 



1 _ ^, 



2(167r2)2 \^47rA/2 J 

t2 rn r — - rj2 

V 1 — T ^ 

24 87r 167r2 



(327r2)2 2e 
(27-f)ln('^^ 



r2 



Tf7 



167r2e 1 24 Stt ^ 



172 
327r2 



In 



/47rT2 



V47rAf2 ) 



1.9785 



, 3A2 
167r2 



V Af2 
V47rAf2 y 



7 - 2r2 



47rry 



A2r2 

1287r2 



5.0669 



ln(l - r2) - ^(2 In 2 - 1) - ^ In^ 



1 



In 2 



2-7 



(30) 



The divergent parts of these diagrams can all be eliminated by a suitable choice of the countertcrm in Eq. (|T^). As 
mentioned above, we use the minimal subtraction prescription. 

III. RESULTS 

Having obtained an expression for the renormalized thermal mass, we can set (5 = 1 and obtain numerical results 
for this mass, rn^ (fj), where fj is determined via the PMS condition: 



drj 



= 0. 



(31) 



ri=V 



At 0{5), r] does not depend on the coupling and so does not generate non-perturbative information, but non- 
perturbative behaviour emerges at 0{S'^). Fig. 3 is a typical plot, showing a clear maximum in rn^ (rf). 
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T=12.50 /Li-0,50 A-0.20 m2 = -l,0 



FIG. 3. Dependence of the thermal mass squared on the variational parameter rj. All masses measured in units of M. 



It will be of particular interest to study the /i dependence of the critical temperature Tc, the signal for the phase 
transition being taken to be ^ = , which leads to an implicit equation for Tc = Tc(^) (with dependence on the 
renormalized bare mass and coupling of the theory suppressed in the notation) . Fig. 4 shows a plot of the contours of 
y for = — M^, A = 0.5 in a region of the (T, ^) plane, with the thicker line indicating the critical temperature. 



m: 



Below this line the LDE breaks down - the thermal mass being a monotonically decreasing function of rj and thus 
lacking extrema. 




1.5 2 



FIG. 4. Contours of the thermal mass squared in intervals of iM in the {T,n) plane for A = 0.5 with m 
axes are in units of M. 



-M\ Both 



In Fig. 5 we present a comparison of Tc{fi) as calculated by the LDE with the first-order estimate 



12 



(32) 



provided by perturbation theory in the high temperature limit 



I' 



FIG. 5. Dependence of the critical temperature on /i at 0(5^) for ni? — 
dotted line indicates the first-order approximation given by Eq. 



-M , A = 0.5, with both axes in units of M. The 



The two curves are quite similar in shape, and converge at large T^. The reason for this behaviour is that in the 



high T limit, it so happens that rj 



which is equivalent to 



1. Examining those contributions from 



114 in ( |29| ) we see that if this is the case, we require the coefficients of the 1/^(1 — r^) terms to disappear, which is 
exactly the same condition as (|3^). 

Fig. 6 illustrates the high-T behaviour of the second-order approximation to the thermal mass, with /i/Tc approach- 
ing a constant as Tc increases; the constant being determined from (^2|) to be \/(A/12). 




FIG. 6. The variation of m/Tc with Tc for = —M^, A = 0.5. The dashed line indicates the constant value -/(A/12). Tc is 
in units of M. 



IV. CONCLUSIONS 



In this paper we have shown how the problem of finite chemical potential for a charged scalar theory can be 
formulated in the context of the linear delta expansion. The graphs to be evaluated are essentially those of ordinary 
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perturbation theory, with modified mass and coupling parameters. The chemical potential appears explicitly in the 
Lagrangian, and also in the Bose-Einstein factors occurring in the momentum integrals. Renormalization has been 
implemented on the lines of Ref . ^ , where the importance of renormalizing before applying the variational aspect of 
the method was emphasized. We obtain unambiguous PMS points, as illustrated in Fig. 3. 

The results have been plotted in Fig. 4 as a contour plot of the thermal mass in the (T, /i) plane, and in Fig. 5 as 
Tc versus fi for given values of m^, A. This latter curve approaches the result of resummed perturbation theory at 
high temperature, but differs significantly from it at lower temperatures. The reason for the convergence at higher 
temperatures can be understood in terms of the properties of the stationary points in the variational parameter ?/. 

As emphasized in the introduction, the problem of a non-zero chemical potential is not amenable to treatment by 
the usual Monte-Carlo lattice method. A lattice version of the present calculation is in progress. 

An extension of the present work which we intend to pursue in the near future is to free it from the dependence 
on the high temperature expansion, thereby enabling us to consider a system with non-zero chemical potential at low 
temperature. 
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APPENDIX A: 




FIG. 7. The two-loop momentum-dependent diagram 



In this appendix we will derive in some detail the expression for the two-loop setting sun diagram. We begin 
by mimicking the approach of Parwani Q to split up the diagram into three parts; containing zero, one and two 
Bose-Einstein factors respectively. Explicitly, 



nr = 



52 



(27r)d-i (27r)'^-i ^^ ^'^^ ' ^^^'^ ^^^'^ ^^^'"^ ' 



(Al) 



where d = 4 — 2e is the space-time dimension and r = p — k — q. If we now use the mixed representation of the 
Matsubara propagator defined in , the frequency sums become trivial and eventually we have 



Hi (zu;„,p) 



<52A2 



(Go (iw„, p) -I- Gi(iw„, p) -I- G2{iuJn, p) ) , 



where 



GoiiuJn,p) = J d[k,q\Sfj,{uJk,UJq,UJr): 



(A2) 



(A3) 



(A4) 
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G2{iuJn,p) = d[k,q] 



+'^k'^q (^Sf,{LOk,UJg,LUr) + {-UJk, LOq, UJr) + S'+(tJfe, -UJq,UJr) - Sf,{LUk,UJq, ~LUr)j 
+nkn^ (^Sfj,{u!k,UJq,iUr) + S'^{-UJk,LUq,U!r) + S^{uJk, -UJq,UJr) ~ Sfj^^LUk, UJq, ~LUr)j 
+ nk'^q (^Sili^k,^q,^r) + S f,{-LUk , UJq, UJr) + S^{u!k, "UJq.LOr) - S^{uJk,L^q, -t^r)) , 



(A5) 



with the definitions 



±{iLUn + /i) + UJk + ^^q + l^r 



Sf,{uJk,^q,^r) = S^{uJk,^q,i^r) + {uJk, UJq , UJr) , 



d [k, q] = M 



4^ d''-ik d-^-iq 



(27r)'^-i (27r)''-i Sw^WgW^ 



We choose to evaluate the self-energy on shell (p — 0,4^,1 = — /^). The reason we evaluate at this point rather 
than at p = Q,iijJn — ^ '\s best explained as follows p6|. In setting up the theory, we are dealing with an effective 
Hamiltonian, HcS = H — fJ.Q rather than the real Hamiltonian, H. We would naturally choose to evaluate the self- 
energy at iujn — f2 for the real Hamiltonian, which means evaluating at icOn = — /i in the effective theory. Following 
this prescription, we find that 



3?[Go(r!-Ai,0)] = - 



2(167r2)2 



1 3-27 2 , 

+ - - In 



(327r 



2\2 



+ (27-¥)ln 



(167r2)2 

We now calculate the contribution from Gi . Writing 

Gi (ia;„ , p) ^G^{iujn,p) + G^{iu!n,p), 





+ 1.9785 


V47rA/2 J 





(A6) 



(A7) 



where G^ represents that part of Gi with an associated BE factor of respectively, and decomposing this factor 
into a UV divergent part and a UV finite part in the manner of the fi = case one obtains 

n[Gt{n - 0)] = ^^0^ + Ft + Ft{n^). 

Here 



^0 — 



and 



with 



F,^{n^) = 



3 


r2 


327r2 


^2 


3 


r2 


327r2 


7r2 


1 





h3{y,±r)-, 
e 

h3{y,±r) 



In 



/47rM2 



+ 2-7 



4(27r) 



dk — - 



dq 

Jo '^q 



q In 



- 6k 



X^ = [il^ - {uJk +t^q +UJk±q)^W^ - (-Wfc +UJq + LJfc±g)^] [(fi ± tJfe)^ - {ujq + UJk±q)'^], 



(A8) 

(A9) 
(AlO) 

(All) 
(A12) 



where the subscript ± refers to the LUk±q term. Following a similar procedure to that of Ref. we find that in the 
high-T limit 



F±(172) 



r2 



1287r2 



In 



0.54597 



(A13) 



Thus, collecting all the terms together. 
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{fl - II, 0)] =Fo + Fi+ ^2(0^), 



(A14) 



with 



3J'2 ^ 



Fi = 



167r4 
167r4 



ht{y,r) 



and 



647r2 



0.54597 



(A15) 
(A16) 

(A17) 



Finally, we consider the contribution from G2, which contains a BE factor for each loop and so is UV finite. We 
write 



with 



4(27r)4 7o cuk Jo 



Y 



±± 



Y_ 



±± 



where 



Y++ = [n^ 



{uJk +UJq + Wfc±,)^][f7^ - (cjfe +U)q- Wfeig)^] [(fi + Wfe - CJ,)^ 
(tJfe +UJq + UJk±q)'^][fl^ - {uJk +(^q- UJk±q)'^][{^ " Wfc + tj,)^ - UJ^ 



- ^k±q\ 



k±qh 
2 , ,2 



(-Wfc + + UJk±q) ] - (Wfe -U>q+ Wfe±g) ] [(1^ - Wfc - U>q) - 



(A18) 



(A19) 



(A20) 
(A21) 
(A22) 
(A23) 



Each of the pieces has a logarithmic IR divergence as J7 — !■ 0. To deal with this, we extract the leading order 

behaviour of In |y_^^/Y^^| in this limit, separating into three terms as follows: 



In 



Y_ 



+± 



In 



~ ±ln2±ln 



kq 



+ |ln 



(A24) 



Collecting those terms involving the In 2 piece alone, which we can take outside the integral, we calculate the coefficient 
of this term in the high-temperature limit a = (30. <C 1, which we call C\n2i to be 



r2\2 



1 + r 



The remainder of K[G2(f^ — 0)] splits into two pieces which we define by 

WV02A 1 r^dk + [^qdq + 



kq 



'kdk 



8(27r)4 7o wfc 



— (n+ +njln 



(A25) 

(A26) 
(A27) 



Extracting the high-temperature limit of if^(n^) gives, after using the symmetry of the integrand under k ^ q to 
restrict the range of the q integration and making the changes of variable fi^ = Wfc, r2?7 = w,. 



o2 /-oo ri 



(A28) 
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where ri^ = l/{e°-^^^^'^ — 1). After making the successive approximations]] 



— rt^ « —2ra 



(e< - 1)2 



(A29) 



in the integrand, we have, for a ^ 1, 



7^2 2 



(A30) 



The calculation of H {Vt ) parallels the case for = 0, leading to 



ij2(n2^ 



3r2 

647r2 



1.50699 



(A31) 
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